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Abstract. I review aspects of the theory of long-duration gamma-ray burst (GRB) central engines. 
I focus on requirements of any model; these include the angular momentum of the progenitor, the 
power, Lorentz factor, asymmetry, and duration of the flow, and both the association and the non- 
association with bright supernovae. I compare and contrast the collapsar and millisecond proto- 
magnetar models in light of these requirements. The ability of the latter model to produce a flow 



with Lorentz factor ^100 while simultaneously maintaining a kinetic luminosity of ~ 10 ergs s _1 
Q • for a timescale of ~ 10 — 100 s is emphasized. 
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^ ; INTRODUCTION 

Models for the central engine of long-duration gamma ray bursts (GRBs) are highly 
constrained by the character of the prompt emission and the afterglow, and — at least in 
^ \ some cases — the fact of an associated supernova (SN). There is little diversity among 

models for the central engine and essentially all can be simply classed as a rotating 
compact object that drives an asymmetric relativistic outflow. 

The "collapsar" model as described in 01 and developed in proposes that GRBs 
arise from the collapse of rapidly rotating Type-Ibc progenitors. A black hole forms 
with an accompanying accretion disk and drives a collimated relativistic outflow along 
the axis of rotation via either neutrino heating or magnetic stresses. 

The "millisecond proto-magnetar" model posits a newly formed rotating neutron star 
(spin period P ~ 1 ms) with surface magnetic field of magnetar strength (B ~ 10 15 G), 
cooling via neutrino radiation on the Kelvin-Helmholtz timescale, ?kh ~ 10 — 100 s, 
and driving a neutrino-heated magneto-centrifugal wind |4|]. Proto-magnetars might be 
produced by rotating Type-Ibc progenitors, the accretion-induced collapse of a white 
dwarf, the merger of two white dwarfs, and/or (potentially) the merger of two neutron 
stars [5]. Thus, they may trace both young and old stellar populations. See also flSS]. 

Here, I compare and contrast some of the basic elements of and requirements on the 
collapsar and proto-magnetar models. I focus on the fact that a viable central engine 
must simultaneously realize a collimated flow with Lorentz factor in the range 100 ^ 
7 ^ 1000, on a 10 — 100 second timescale, while producing a kinetic luminosity of 
~ 10 50 ergs s . While many models have been proposed that can in principle meet 
these requirements, the proto-magnetar model realizes them quantitatively. 



THE PROGENITOR: ANGULAR MOMENTUM 



Collapsar. The model is that of a central black hole, formed by core-collapse of 
a massive star, with a centrifugally-supported disk. The latter requirement sets the 
minimum angular momentum required for the collapsar mechanism: y' m i n ~ 1.5 x 
10 16 M3 cm 2 s _1 is required for disk formation at the ISCO of a maximally -rotating 
M3 = M/3 M & black hole. For a disk that extends to larger radius and/or a non-maximal 
black hole, j m i n increases. The models of [@] that produce a strong "Ni wind," which is 
needed to power the bright associated Type-Ibc SN lightcurve have j ~ 10 17 cm 2 s _1 . 

Proto-Magnetar. The proto-magnetar model does not require disk formation. In- 
stead, it posits the existence of a neutron star with a large-scale magnetic field of strength 
B ~ 10 15 G, mass Mi .4 = M/1.4M Q , radius R[q = i?/10km and spin period of order 
P\ = P/l ms 0]. In a detailed set of calculations, [@] find that (depending on B) P ^ 2 ms 
is required to produce conditions favorable for a GRB, impling a minimum specific an- 
gular momentum for the pre-collapse iron core of j m i n ~ 3 x 10 15 R 2 Q P^ 1 cm 2 s _1 . This 
is a factor of 5 smaller than j' m i n for the collapsar model, and a factor of ~ 30 smaller than 
the model of [2] that launches a disk wind. In absolute terms, a spin frequency (Q.) cor- 
reponding to P ~ 2 ms is relatively small, just ~ 25% of breakup, ^ max ~ (GM/7? 3 ) 1 / 2 . 

The additional requirement on the proto-magnetar model is that the neutron star must 
generate or be born with a magnetar- strength field. That large B might accompany small 
P was proposed by [floLITlll . who argued that when P becomes shorter than the convective 
overturn timescale (Rossby number < 1) neutron stars undergo strong dynamo action, 
producing a magnetar. A weakness of the proto-magnetar model is that a strong poloidal 
field is assumed to exist without being generated self-consistently (e.g., [|4j, li2D). In 
absolute terms, the energy associated with a 10 15 G field is much less than the rotational 

energy (5 eq ,rot ~ [AitpR 1 ^ 2 ] 1 ' 2 ~ 2 x lO^pJf * lo^f 1 G, where p 14 = p/10 14 g cm" 3 ) 
of the neutron star and/or the total energy carried in convective motions during ?kh- 

Summary. j m [ n for collapsars is significantly larger than for proto-magnetars. If a 
progenitor stellar population evolves in such a way as to produce conditions favorable 
for a collapsar, it is hard to see how it should not produce progenitors with smaller 
specific angular momentum in the range needed to power proto-magnetar-driven GRBs. 
Perhaps in this way collapsars represent an externum of the GRB population and proto- 
magnetars more continuously connect with the SN population. 

THE FLOW: KINETIC POWER & LORENTZ FACTOR 

In the internal shock model, it is important that the flow that generates the GRB achieve 
100 1000 while simultaneously producing a kinetic luminosity of ~ 10 50 ergs s _1 

on a t ~ 10— 100 s timescale. Reference \ j l3\ \ writes that "A nagging question in all 
these models is what produces the 'observed' ultra-relativistic flow? How are 10 5 M 
of baryons accelerated to an ultra-relativistic velocity with 7 ~ 100 or larger? Why is 
the baryonic load so low? Why isn't it lower? There is no simple model for that. An 
ingenious theoretical idea is clearly needed here." See the lucid and concise discussion 
in [ | 1 4T| of the physics that sets the upper and lower limits on 7. 

Collapsar. There have been few attempts to derive the Lorentz factor self- 
consistently in the collapsar model. Reference |2J] showed that vv annihilation above 



the poles of the black hole results in efficient energy deposition, generating very large 
specific entropy and high relativistic enthalpy. It is natural that the flow should then 
accelerate to relativistic velocities, as in II 151 LL6Q. A collimated relativistic flow might 
also be accelerated via magnetic stresses along the rotational axis, producing a highly 
Poynting-flux-dominated jet. Calculations of the global structure of accretion disks 
indicate that the polar funnel does realize very large magnetization, a, the magnitude of 
which is limited only by numerics as the pole becomes essentially baryon-free II 1711 . It 
has been suggested that cross-field neutron diffusion might set the asymptotic Lorentz 
factor in an otherwise baryon-free jet lfl8|1 : however, see Il9ll . 

There is little doubt that the flow above the pole of a rotating black hole will become 
highly relativistic as a result of either neutrino energy deposition or magnetic stresses, 
or both. The question is, does the collapsar model simultaneously produce a jet with 
100 ^ 7^ 1000 and E ~ 10 50 ergs s _1 ? What is the time-evolution of y, and how 
does it relate to the time-dependence of the accretion rate? These questions are as yet 
unanswered and they represent the largest single gap in the collapsar model of GRBs. 

Proto-Magnetar. Regardless of the mechanism of core-collapse SNe, neutron stars 
are born hot — with a central temperature of 10's of MeV — and they radiate their 
gravitational binding energy in neutrinos as they cool, contract, and deleptonize on the 



Kelvin-Helmholtz timescale of ?kh ~ 10— 100 s [20, |21|]. During this epoch, neutron 
stars drive thermal neutrino-heated hydrodynamical winds (lA l23|l . For a neutron star 
with M = 1.4M©, R = 10 km, rotation frequency Q., and and a magnetar- strength mag- 
netic field, the mass-loss rate during ?kh is given approximately by ll23l.l4l.l9ll 

M(t) « 10- 6 4( 2 51 (0 exp[Q(0 2 /^c] M s- 1 , (1) 

where the total neutrino luminosity (L v ) is indexed by the v e luminosity Ly e ,5i = 
Ly e /10 51 ergs s _1 (typically L v ~ 5Ly e ) and it has been assumed that Ly e oc (£v e ) 4 , 
where (£y e ) is the average energy. The normalization of eq. (OQ) and its L v -dependence 
follow from the physics of the weak interaction, and, in particular, the cross-section 



for the charged-current processes v e n «-> pe~ and v e p <-> ne + 11221 12311 . This physics 
is a significant part of the answer to the questions posed by jl3|l at the beginning 
of this section in the proto-magnetar model. The exponential factor in equation (OQ) 
accounts for the enhancement by magneto-centrifugal forces when B and Q. are large 
(£2 C « 2300L?'°g 1 rad s 1 ; ^\). For P ^ 2 ms this factor becomes important 0]. 

One of the most important components of the prot o-magnetar model is that as the 
neutron star cools L v decreases on a timescale ?kh I2JJ,LZL|4|]. As a result of equation 
(OQ), M decreases concomitantly. For this reason, for fixed surface magnetic field strength, 
the wind becomes increasingly magnetically-dominated and relativistic as a function of 
time. This is quantified by the magnetization at the light cylinder (R^ = c/Q. ~ 50Pi km): 

o L (t) =B{R L ) 2 /[ATtp(R L )c 2 ] ocM(ty l . (2) 

At early times of order ~ 1 s after the SN shock is launched L v is high enough that the 
wind is primarily driven by neutrino heating, M is large, Ol is less than unity, and the 
flow is non-relativistic. As L v decreases, the wind becomes increasingly magnetically- 
dominated and it transitions to non-relativistic, but magneto-centrifugally dominated. On 



a few- second timescale, Ol becomes larger than unity and the wind becomes Poynting- 
fiux dominated. It becomes increasingly so on a timescale ?kh- If we assume that there 
is efficient conversion of the magnetic energy to bulk kinetic energy (via, e.g., magnetic 
dissipation; see [24]), then the asymptotic Lorentz factor of the wind is y ~ Ol- 



As an example, in the models of [9] (see, e.g., Fig. 1 of 112511 ). we find that Oi ~ 100 at 
t ~ 20 s, wheni? « 10 50 ergs s , for a proto-magnetar with 5 = 3 x 10 15 G and/ 5 = 1 ms. 
Thus, proto-magnetars drive relativistic flows with the kinetic luminosity, timescale, and 
Lorentz factor appropriate for producing GRBs. 

The time evolution of the system is governed by perhaps the single most significant 
input to the proto-magnetar model: the time evolution of L v and (e v ) : as L v decreases, M 
decreases (eq. CO, Ol increases (eq. [2]), and, because E scales with a postive power of M 
ACS 111], it decreases. The cooling epoch eventually ends, L v and M drop precipitously, 



and the proto-magnetar transitions to more "pulsar"-like. As recently shown by lll9h . the 
fact that <7l(0 ~ 7(0 increases monotonically with time throughout ?kh implies high 
radiative efficiency in the internal shock model for the prompt emission. 

In J3, & Qji] we have taken L v (t) and (£ v )(t) from ll2lll . informed by the early- 



time calculations of, e.g., I126ll . However, the models of [2lj] are for non-rotating non- 
magnetic proto-neutron stars and it is possible that models including these effects (more 
appropriate for proto-magnetars) will change L v (t) and (£ v )(t) quantitatively. 

Summary. The ability of the proto-magnetar model to link L v , M, E, and Ol ~ J 
and their mutual time-dependence is a significant strength of the model and may allow it 
to be ruled out or (potentially) confirmed for some subsets of the global GRB population. 

THE EXPLOSION: (NON-)ASSOCIATION WITH SUPERNOVAE 



Several events link GRBs and SNe directly (e.g., GRB 030329 & SN 2003dh; [27]). 
Conversely, there are cases for which firm limits on an associated SN (M[ 56 Ni] < 
10~ 3 - 1O~ 4 M ) can be established (e.g., 060505, 060614; SSH). 

Collapsar. The collapsar model provides an explanation for both the association and 
non-association with SNe. High angular momentum cores produce extended collapsar 
disks, which launch non-relativistic winds as in ||2j]. Calculations taking into account the 
electron fraction (Y e ) changing charged-current interactions show that even though the 
midplane of the disk is neutron-rich, the outflow becomes de-neutronized with Y e ps 0.5, 
which then produces the 56 Ni needed to power a SN lightcurve. In contrast, low angular 
momentum cores that do not produce extended disks might not drive Ni winds, although 
a relativistic jet above the poles is potentially possible. Both scenarios require a young 
stellar population. This is qualitatively different from the proto-magnetar model. 

Proto-Magnetar. Some proto-magnetars may be formed from the collapse of a 
Type-Ibc progenitor. In this scenario, one imagines that the SN mechanism in some 
form (e.g., the "neutrino mechanism" rf3lll ) launches a SN shock with ~ 10 51 ergs that 
produces 56 Ni via explosive nucleosynthesis. This would explain the presence of an 
accompanying SN, but not the bright character of some GRB-SNe. If Ni in excess of 



that produced in non-GRB Type-Ibc SNe is indeed required in some cases (see [|32|D 
there are two ways in which it might be produced by proto-magnetars: (1) some of 
the initial rotational energy of the core may be tapped rapidly via magnetic stresses, 
enhancing the SN shock energy as it is launched L26, 13311 and/or (2) as the initial slow 



SN Shockwave is moving outward it is shocked by the subsequent, highly-energetic 
proto-magnetar wind, again enhancing the shock energy Ho. Depending on the angular 
distribution of the wind kinetic energy, the latter option requires that ^ 10 51 ergs is 
extracted from the proto-magnetar on a ^ 1 — 2 s timescale. 

Proto-magnetars (and their GRBs) may also be formed by the accretion-induced 
collapse of white dwarfs and the merger of white dwarfs. In this scenario there is no 
explosive nucleosynthesis, essentially no 56 Ni yield, and no accompanying SN [l9ll. 
Proto-magnetar GRBs could then potentially trace both relatively young and old stellar 
populations. This is qualitatively different from the collapsar model. 

Summary. The 56 Ni yields of non-GRB Type-Ibc SNe are presumably generated 
by explosive nucleosynthesis. It would be remarkable if an entirely different mechanism 
— i.e., the Ni wind from a collapsar disk — was responsible for 56 Ni production in 
GRB-SNe given the quantitative similarlity between the 56 Ni yield distributions of GRB- 
SNe and non-GRB-SNe |[32(l . Yet, this possibility is not excluded. A potentially more 
natural explanation is that the similarity in 56 Ni yields between GRB-SNe and non-GRB- 
SNe evidences the same underlying physical mechanism. The proto-magnetar model 
provides a simple and natural way of understanding this, as well as the association and 
the non-association of SNe with some GRBs. It further allows for a simple interpretation 
of the continuum in properties between the explosive events of non-GRB-SNe and the 
GRB population itself (see, e.g., Fig. 5 of Q34fl ): the diversity in B and Q. for neutron stars 
at birth leads to a diversity in the amount of matter accelerated to relativistic velocities 
and the energy and distribution of energy as a function of 7 for the explosion as a whole. 

ASYMMETRY: BEAMING & COLLIMATION 

Collapsar. Because of the geometry of the system, the collapsar model provides 
a natural explanation for how the flow becomes highly-collimated, whether driven by 
neutrino energy deposition J2J] or magnetic stresses |l7|1 . 

Proto-Magnetar. The kinetic luminosity of a relativistic Poynting-flux-dominated 
proto-magnetar wind emerging into vacuum is distributed broadly around the equatorial 



plane 111211 : production of a jet is not trivial. However, work by 11251 . 13511 shows that 
the interaction between the outflow and both the overlying stellar progenitor and the 
preceding SN shock acts to tightly collimate the relativistic flow. The mechanism was 
suggested in the GRB context by ll3~6ll . based on models of [H3]. At radii much larger 
than Rl, the toroidal component of the wind magnetic field (B^) is much larger than its 
polodal component. The wind shocks on the exploding stellar envelope and produces a 
relativistically hot, quasi-hydrostatic bubble. If B§ is large enough, the bubble expands 
in the polar direction as a result of both hoop stress at the equator and the confinement of 
the bubble by the overlying stellar envelope. The result is a relativistic jet. The models 
of 113 511 indicate that little of the incident wind energy is coupled to the "spherical" 
component of the explosion (the SN) and, therefore, the asymptotic Lorentz factor of 
the jet reflects the Lorentz factor of the incident proto-magnetar wind. 

Summary. Based on analogy with observed jets in AGN and X-ray binaries it may 
be argued that the picture of collimation in the collapsar scenario is more natural than 



that for the proto-magnetar. Nevertheless, [25, 35] provide a compelling picture of proto 



magnetar wind collimation. A major focus of future work on the proto-magnetar model 



will be to understand the interaction of the wind on the progenitor and the potential 
feedback on the wind itself if, for some parameters, the reverse shock propagates inside 
the fast, Alfven, and slow magnetosonic surfaces H35I1 . 
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